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In this paper we perform a global analysis of the constraints on the inflationary parameters in the 
presence of dynamical dark energy models from the current observations, including the three- year 
Wilkinson Microwave Anisotropy Probe (WMAP3) data, Boomerang-2K2, CBI, VSA, ACBAR, 
SDSS LRG, 2dFGRS and ESSENCE (192 sample). We use the analytic description of the infla- 
tionary power spectra in terms of the Horizon- flow parameters {ti}. With the first order approxi- 
mation in the slow-roll expansion, we find that the constraints on the Horizon-flow parameters are 
ei < 0.014 (95% C.L.) and e 2 = 0.034 ± 0.024 (la) in the ACDM model. In the framework of 
dynamical dark energy models, the constraints become obviously weak, e\ < 0.022 (95% C.L.) and 
€2 — —0.006 ± 0.039 (lcr), and the inflation models with a "blue" tilt, which are excluded about 
2a in the ACDM model, are allowed now. With the second order approximation, the constraints 
on the Horizon-flow parameters are significantly relaxed further. If considering the non-zero 63, the 
large running of the scalar spectral index is found for the ACDM model, as well as the dynamical 
£^ ■ dark energy models. 
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\0 • I. INTRODUCTION 

Inflation in the very early universe is the most attractive paradigm, which is driven by a potential energy of a 
43 ■ scalar field called inflaton and its quantum fluctuations turn out to be the primordial density fluctuations which seed 
Q-r the observed large scale structures (LSS) and the anisotropics of cosmic microwave background radiation (CMB). 

Inflation theory has successfully passed several non- trivial tests. The current cosmological observations are in good 
?-h ' agreement with an adiabatic and scale invariant primordial spectrum, which is consistent with single field slow-roll 
inflation predictions. And the large angle anti-correlation is found in the temperature-polarization power spectrum, 
which is the signature of adiabatic superhorizon fluctuations at the time of decoupling [l| . 

In 1998, the analysis of the redshift— distance relation of type la supernova (SNIa) revealed the existence of the 
£SJ ■ another stage of accelerated expansion that started rather recently when a mysterious new energy component dubbed 
' dark energy (DE) dominated the energy density of the Universe [2| . The nature of dark energy is among the biggest 
problems in modern physics and has been studied widely. The simplest candidate of dark energy is the cosmological 
constant (CC) however it suffers from the fine-tuning and coincidence problems Q. To ameliorate these dilemmas 
some dynamical dark energy models such as Quintessence [H , Phantom jf| and K-essence @ . Given our ignorance of 
the nature of dark energy, constraining the evolution of DE the equation of state (EoS) by cosmological observations 
is of great significance. Interestingly, there exists some hints that the EoS of dark energy has crossed over —1 at least 
once from current astronomical observations @, namely Quintom dark energy model, which greatly challenges the 
above mentioned dark energy models. 

In 2006, the WMAP group @ obtained the constraint on the scaler spectral index n s = 0.958 ± 0.016, which 
deviates from the simple scale-invariant primordial spectrum and disfavors the inflationary models with a "blue" tilt 
at more than 2a. Alternatively, the scale- invariant Harrison-Zel'dovich-Peebles (HZ) spectrum (n s = l,r = 0) is 
disfavored about 3a Q . And the large running of the scalar spectral index is still allowed [l(| ■ It seems that the 
scale-invariant spectrum is disfavored and the dynamics of Inflation has been detected. Similar results have also been 
found in the literature from the current observational data [111 ]. But it's noteworthy that these analysis are based on 
the ACDM model. In the framework of dynamical dark energy models, the constraints on the inflationary parameters 
can be relaxed due to the degeneracy among the inflation and dark energy parameters [eH [l3| . 

In this paper we use the current cosmological observations to carry out a first detailed study on the inflationary 
parameters in terms of the Horizon-flow parameters {e^} in the presence of dynamical dark energy models. Our results 
show that the dynamics of dark energy models weaken the constraints on the Horizon-flow parameters significantly. 
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II. METHOD AND DATA 



In order to compare the theoretical predictions of inflation models with the cosmological observations, we often 
parameterize the primordial power spectra of scalar and tensor perturbations as: 



P s (k) = A s exp 
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where rij is the spectral index, a; denotes the running of the spectral index and fc* is the pivot scale. In this paper 
we use the analytic description of the inflationary power spectra in terms of the Horizon- flow parameters {f j} 1 , which 
are based on the Hubble parameter during inflation and its derivatives, defined as [l4l |: 
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where N is the number of e- foldings. Following the Eq.|T]), t he p ower spectrum can be obtained as an expansion of 
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the power spectrum in terms of the logarithmic wavenumber [16j |: 
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where P s q — H 2 G/tt€i, P t o = 16H 2 G/tt and the coefficients bi given in Ref.flq] are related to the Horizon-flow 



parameters {ei} ■ The coefficients for the scalar spectrum are: 

_2 



b s0 = -2(C + l)ei - Ce a 

2 



-2C-7 



III 

12 



- C 2 - 3C - 7 eie 2 



- 1 



7T 

b s i = n s - 1 = — 2ei - e 2 
6 S 2 = a s = -2eie 2 - £2^3 

and those for the tensor spectrum are: 
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where C = ln2 + 7£ — 2« —0.7296 (je is the Euler-Mascheroni constant). And the ratio of amplitudes of the scalar 
to the tensor at the pivot scale is: 
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At the first order approximation Eq. (|10p becomes the well-known consistency relation of inflation r = —8n t . 

For dark energy, we choose the commonly used parametrization of the dark energy equation of state (EoS) as 1£ 



WDE(a) = Wo +wi(l- a) 



(11) 



where a = 1/(1 + z) is the scale factor and w\ = —dw/da characterizes the "running" of the equation of state (RunW 
henceforth) . For comparison we also consider the ACDM model and the dark energy model with a constant equation of 
state (WCDM henceforth). When using the MCMC global fitting strategy to constrain the cosmological parameters, it 
is crucial to include dark energy perturbation [9l Il9l. |20| . However, it is divergent when the parameterized EoS crosses 
w = — 1 (2lj . By virtue of Quintom dark energy model Q, whose EoS can smoothly cross w = — 1, the perturbation 
at the crossing points is continuous. Thus we have proposed a technique to treat dark energy perturbation in the 
whole parameter space, including w > — 1, w < — 1 and at the crossing points. For details of this method, we refer 
the readers to our previous companion papers (20j . 



1 In the literature [l5ll the slow-roll parameters, e, rj and £, are also used to constrain the inflation models. 

2 Ref. Il7ll used the method of comparison equations in the study of the cosmological perturbations and obtained the similar coefficients 
bi. 
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TABLE I. Mean lcr constraints on the cosmological parameters using the current observations. For columns I and II, the Horizon-flow 
parameters are obtained at the first and second order approximation in the slow-roll expansion respectively. For the weakly constrained 

parameters, ei and r, we quote the 95% upper limit instead. 
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In this study, we have modified the publicly available Markov Chain Monte Carlo package CAMB 3 / CosmoMC 4 
[23l ] to include the dark energy perturbation and the public available code by Leach and Liddle for the primordial 
spectrum 5 [24|. We assume purely adiabatic initial conditions and a flat universe. Our most general parameter space 
is: 

V = u c , Q s , t, w , wi,e 1 ,e 2 ,e 2 e 3 , log[10 10 A s ]) , (12) 

where u>b = Qbh 2 and u> c = £l c h 2 are the physical baryon and cold dark matter densities relative to the critical density, 
S is the ratio (multiplied by 100) of the sound horizon to the angular diameter distance at decoupling, r is the optical 
depth to reionization, A s is defined as the amplitude of the primordial spectrum. For the pivot scale of the primordial 
spectrum we set A* = 0.05 Mpc -1 . 

In the computation of CMB we have included the WMAP3 Temperature- Temperature (TT) and Temperature- 
Polarization (TE) power spectra with the routine for computing the likelihood supplied by the WMAP team 0, [2f| 
as well as the smaller scale experiments, including Boomerang-2K2 (26|, CBI [13], VSA [28[ and ACBAR [29[. For the 
Large Scale Structure information, we have used the Sloan Digital Sky Survey (SDSS) luminous red galaxy (LRG) 
sample [3(| and 2dFGRS [3l|. To be conservative but more robust, in the fitting to the SDSS LRG sample we have 
used the first 15 bins only, 0.0120 < /c cu < 0.0998, which are supposed to be well within the linear regime. For SNIa 
we have marginalized over the nuisance parameter & The supernova data we use are the ESSENCE (192 sample) 
data [33|]. Furthermore, we make use of the Hubble Space Telescope (HST) measurement of the Hubble parameter 
Hq = lOOh km s _1 Mpc -1 [HI by multiplying the likelihood by a Gaussian likelihood function centered around 
h = 0.72 and with a standard deviation a = 0.08. We also impose a weak Gaussian prior on the baryon density 
fift/i 2 = 0.022 ± 0.002 (1 a) from Big Bang Nucleosynthesis (35|. Simultaneously we will also use a cosmic age tophat 
prior as 10 Gyr < to < 20 Gyr. 

For each regular calculation, we run 8 independent chains comprising of 150,000 — 300,000 chain elements. The 
average acceptance rate is about 30%. We test the convergence of the chains by Gelman and Rubin criteria [36[ and 
find R — 1 is of order 0.01 which is more conservative than the recommended value R — 1 < 0.1. 



III. RESULTS 

We summarize our main global fitting results in Table I. Table I lists all of the relevant one-dimensional median 
values and la constraints. Shown together are the corresponding reduction of Xmin values compared with the ACDM 
model. For the constraints on ei and r only 2a upper bounds have been shown. 

Firstly we consider the first order approximation in the slow-roll expansion, where the relevant parameters are ei 
and £2- In this case the running of the spectral index vanishes. In the ACDM model, illustrated in the left up panel 



3 http://camb.info/. 

4 http:/ /cosmologist.info/cosmomc/. 

5 http: / /astronomy. Sussex. ac.uk/ ~sleach/inflation/camb_ inflation. html/. 
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FIG. 1: Constraints on the inflationary parameters at first order in the slow-roll approximation in the ACDM (blue dash-dotted 
lines), WCDM (red dash lines) and RunW (black solid lines) dark energy models respectively. The left up panel and right 
down panel are the one dimensional marginalized distribution of Horizon-flow parameters ei and £2 ■ The left down panel gives 
the two dimensional constraint on (e2,ei). And the right up panel gives the two dimensional constraint on (n a ,r). The contours 
stand for the 68% and 95% confidence level. The two solid magenta lines delimit the three classes of inflation models, namely, 
small-field, large-field and hybrid models. The blue points are predicted by m 2 (j) 2 model and Xcj> 4 model respectively with the 
number of e-foldings, N, being 50 — 60 for m 2 (fi 2 model and 64 for \<f> 4 model. 



and right down panel of Fig[TJ we obtain the constraints on the Horizon-flow parameters are t\ < 0.014 (95% C.L.) and 
e 2 = 0.034±0.024 (la) and consequently with Eq.© and Eq.QI)]) we obtain the spectral index n s = 0.955±0.019 (la) 
and the tensor-to-scalar ratio r < 0.223 (95% C.L.), which is in good agreement with the WMAP's results @. 

From the right up panel of Fig[TJ one can see that a pure HZ spectrum for scalar perturbations with no tensors 
(n s — 1, r = 0) is clearly disfavored at more than 2a by the current observations. We plot the same constraints in 
terms of the Horizon- flow parameters ei and e 2 in the left down panel of Fig[Tp. For the simple monomial chaotic 
models, the single slow-rolling scalar field with potential V((f>) ~ m 2 </> 2 , which predicts (n S) r) = (1 — 2/N,8/N), 
is well within la region, while another single slow-rolling scalar field with potential V(<t>) ~ X<f> , which predicts 
(n s ,r) = (1 - 3/N, 12/JV), is excluded by more than 2a in the ACDM model [H H [H . 

However, the current observational data don't exclude the dynamical dark energy models and especially mildly 
favor a class of models with EoS across the cosmological constant boundary Due to the degeneracy between 
inflation and dark energy, it's necessary to perform an analysis of global fitting allowing simultaneously the dynamics 
in both inflation and the dark energy sector. 

In the framework of dynamical dark energy models, we find that the constraints on the Horizon- flow parameters 
ei and e 2 and the derived parameters n s and r have been weakened dramatically as shown in the one dimensional 
distribution plots of Fig[T] Quantitatively, the constraints on the Horizon-flow parameters are t\ < 0.022 (95% C.L.) 
for both models and e 2 = 0.004 ± 0.036 (la) for WCDM model, e 2 = -0.006 ± 0.039 (la) for RunW model. And 
constraints on the derived parameters are n s = 0.977 ± 0.024 (la) for WCDM model, n s = 0.986 ± 0.028 (la) for 
RunW model, and r < 0.35 (95% C.L.) for both models. The mean value of n s gets closer to n s — 1 and the 95% 



Ref. [37| constrained on the dynamics of Inflation in the (e2,ei) plane straightforwardly. 



5 



upper limit of r is relaxed due to the degeneracy that the tensor fluctuation and the dark e nergy co mponent, through 
the ISW effect, mostly affect the large scale (low multipoles) TT power spectrum of CMB [12l. l44j] . 

Because of this degeneracy, in the two dimensional plot of FigJU we can find that the allowed parameter space 
is enlarged dramatically. Consequently the HZ spectrum, disfavored about 3a in the ACDM model, can be allowed 
within the 2a region in the presence of the dynamics of dark energy. And interestingly many hybrid inflation models, 
excluded in the ACDM model, revive in the framework of dynamical dark energy models as illustrated in FigQ] 13]. 




0.02 0.04 -0.1 0.1 0.2 0.1 0.2 




0.85 0.9 0.95 1 1.05 -0.2 -0.1 0.2 0.4 0.6 



n a r 
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FIG. 2: Marginalized posterior probability distributions for the Horizon-flow parameters, ei, £2 and £2^3, and the derived 
inflationary parameters, n a , a 3 and r, up to second order approximation in the slow-roll expansion in the ACDM (blue dash- 
dotted lines), WCDM (red dash lines) and RunW (black solid lines) dark energy models respectively. 



With the second order approximation in the slow-roll expansion, one has to consider the third Horizon-flow param- 
eter 63. In our calculation we choose 62^3 as the basic parameter directly instead of e 3 . Practically this will make our 
numerical calculation much more efficient. And we notice that the second order formalism are valid in the limit of 
£2^3 <§C 1) but notes <C 1 [H, [45|. We have also checked with £3 as a parameter in the fitting and found the constraint 
is very poor [HiS S|. 

In Fig|2]we show the one dimensional marginalized posterior probability distributions for the Horizon-flow param- 
eters, ei, 62 and 6263, and the derived inflationary parameters, n s , a s and r, up to second order approximation in the 
slow-roll expansion. The constraints on the Horizon-flow parameters become weaken obviously in the dynamical dark 
energy models relative to the ACDM model. 

Another result is the appearance of large running of scalar spectral index which has been found in the literature 
0, [Io| . From the current observations, illustrated in Fig[2J we obtain the fully marginalized value of £263 and a s : 
e 2 e 3 = 0.089 ± 0.042, a s = -0.092 ± 0.045, which deviate from zero with more than 2a, in the ACDM model. This 
large value of running a s violates the inequality: 
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for k far away from the pivot scale fc* . One possible explanation is that current observations are not accurate enough 
to determine the scalar spectral index yet. This large running needs much more observation data, such as PLANCK 



measurement 13j, to verify. Indeed, these constraints may be affected by the Lyman-a forest data [411 ]. However, if 
this large value of running a s can be confirmed by the future measurement, this would be a great challenge to the 
single-field inflation model described by the slow roll expansion which can not produce this large, negative running of 
scalar spectral index a s [42| . At that time, two or more Inflationary history or breaking down the slow roll expansion 
would be needed [43| . 

In the framework of dynamical dark energy models, the large running of scalar spectral index is still allowed. The 
constraints on a s are a s — —0.082 ± 0.045 and a s = —0.075 ± 0.046 for the WCDM and RunW models respectively. 
The mean value of a s slightly shift and the error bars are unchanged. It seems that the correlation between dark 
energy parameters and the running a s is weak [44[ . This weak correlation might be understood from the distinct 
effects on CMB TT power spectrum. For the dark energy parameters Wq and Wi, the effect on CMB TT power 
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spectrum can be somewhat identified with a constant effective equation of state [461 ] : 

Jdafl(a)w(a) 

WcS = Jdailia) ■ (W) 

Keeping other cosmological parameters unchanged, the TT power spectrum will be shifted to larger scalar if the 
effective EoS w e s becomes larger. However, the negative running a s will suppress the amplitude of the spectrum, 
which can not be mimicked by adjusting the dark energy parameters only. This may be the reason of the weak 
correlation between the dark energy parameters and the running of scalar spectral index a s . 



IV. SUMMARY 



In this paper we perform an analysis of global fitting on the inflationary parameters in terms of the Horizon-flow 
parameters {e^} in the presence of dynamical dark energy models from the current observations. Our analysis shows 
that the dynamics of dark energy generally weakens the constraints on inflationary parameters, due to the degeneracy 
between dark energy and inflation parameters. 

With the first order approximation in the slow-roll expansion, the constraints on e± and ti can be significantly relaxed 
and the allowed parameter space of (e2,ei), {n s ,r) panels are enlarged relative to the ACDM model. Consequently 
the HZ spectrum (n s = 1, r = 0), disfavored about 3<r in the ACDM model, can be allowed within 2a in the presence 
of the dynamics of dark energy. Interestingly many hybrid inflation models, especially for models with a "blue" tilt 
(n s > 1), excluded in the ACDM model, revive in the framework of dynamical dark energy models. 

With the second order approximation in the slow-roll expansion, we use the parameter de?/dN = £2^3 to do the 
calculations instead of £3. We find that the constraints on the Horizon-flow parameters become weakened further and 
the large running of scalar spectral index is still allowed, in the framework of dynamical dark energy models. The 
degeneracy between the running of the scalar spectral index and the dynamics of dark energy is weak. 
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